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Abstract 

In a previous paper, we have given an algebraic model to the set of intervals. Here, we apply this 
model in a linear frame. We define a notion of diagonalization of square matrices whose coefficients 
are intervals. But in this case, with respect to the real case, a matrix of order n could have more 
than n eigenvalues (the set of intervals is not factorial). We consider a notion of central eigenvalues 
permits to describe criterium of diagonalization. As application, we define a notion of Exponential 
mapping. 



1 The associative algebra ER 

In p] , we have given a representation of the set of intervals in terms of associative algebra. More precisely, 
we define on the set IR of intervals of R a R- vector space structure. Next we embed IR in a 4-dimensional 
associative algebra. This embedding permits to describe a unique distributive multiplication which 
contains all the possible results of the usual product of intervals and the monotony property is always 
conserved. Moreover, this new product is minimal with respect the distributivity and the monotony 
properties. 

In this section, we present briefly this construction (for more details, see [I]). 

1.1 Vector space structure on IK 

Let IR be the set of intervals of R, that is 

IR = {[a,b],a,be R, a < b}. 

This set is provided with a semi-group structure that we can complete as follow: we consider the 
equivalence relation on IR x IR: 

(X, Y) ~ (Z, T) <^> X + T = Y + Z 

for all X, Y, Z,T e IR. The quotient set is denoted by IR. The addition of intervals is compatible with 
this equivalence relation : 

{X, Y) + (^i) = (x + z,y + t) 

where (x, y) is the equivalence class of (X, Y). The unit is = {(X, X), X e IR} and each element has 
an inverse 



Then (IM, +) is a commutative group. We prove also in [I], that any equivalence class admits a canonical 
representant of type (K , 0) or (0,K) or (a, a), with K £ ffi and a = [a, a], a £ K. We provides the 
group (IM, +) with a real vector space structure, the external product being given by 

f a • (if, 0) = (aif,0) 
{ a ■ (0,K) = (0, aK) 

if a > 0. If a < we put /3 = —a and 

1 a- (0,10 = (AM) 
The triplet (IR, +, •) is a real vector space. 

Remark. To simplify notations, we write (K,0) or (0, K ) in place of (1^,0) or (0, K). 
1.2 The associative algebra A4 

Recall that by an algebra we mean a real vector space with an associative ring structure. Consider the 
4-dimensional associative algebra whose product in a basis {ei, e%, e$, e^} is given by 
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The unit is the vector e 4 + e2- This algebra is a direct sum of two ideals: A4 = I\ + I2 where Ii is 
generated by ei and e 4 and I2 is generated by e2 and 63. It is not an integral domain, that is, we have 
divisors of 0. For example e 4 ■ e% = 0. The ring A4 is principal that is every ideal is generated by one 
element. The cartesian expression of this product is, for x = {x\, X2, £3, X4) and y — (2/1,2/2,2/3,2/4) in 
A4: 

X ■ V = (Xiyi + X4V4, X 2 V2 + £32/3, £32/2 + Z22/3, £42/1 + £13/4). 

The multiplicative group A\ of invertible elements is the set of elements x = (x\, x%, X3, £4) such that 

( x 4 ^ ±Xl, 

\ x 3 ^ ±x 2 - 

If x £ A\ we have: 

-l _ ( X l X 2 £3 Xj \ 

\ ry>2 ^, 2 ' ™ 2 ryi2 ' ™2 /y> 2 ' ™2 T 1 / 

4 *** 2 3 2 3 1 4 / 

1.3 A product in M 

We define a correspondence between IM and A. Let tp be the map 

ip : M — > Aa 
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given by 



(xi , x 2 , 0, 0) if Xi , x 2 > 0, 
y(K, 0)= { (0, x 2 , -Xi, 0) if Xi < and x 2 > 0, 
(0,0, -a;i, -x 2 ) \ix 1 ,x 2 < 0, 



and 

for any if = [xi,x 2 ]. Then we have 



lp(0,K) = -lp(K,0). 



— xi, — x 2 , 0, 0) if xi > 0, 
^(0,if) = <( (0,-x 2 ,-xi,0) if xix 2 < 0, 
(0, 0, -xi, -x 2 ) if x 2 < 0. 

Thus the image of M in A\ is constituted of the elements 

(xi,x 2 ,0, 0) with < xi < x 2 which corresponds to ([xi,x 2 ],0), 
(0,x 2 ,— Xi,0) with x\ < < x 2 which corresponds to ([xi,x 2 ],0), 
(0, 0, — Xi, — x 2 ) with Xi < x 2 < which corresponds to ([xi, x 2 ], 0), 
(— Xi,— x 2 ,0,0)with < xi < x 2 which corresponds to (0, [xi,x 2 ]), 
(0, — x 2 ,xi,0) with x\ < < x 2 which corresponds to (0, [xi,x 2 ]), 
(0,0, Xi,x 2 ) with xi < x 2 < which corresponds to (0, [xi,x 2 ]). 



x\-V\= x 3 - 2/3, 
x 2 - 2/2 =x i -y i 



The map ip : IR — >Ai is not linear. We introduce in A4 the following equivalence relation 1Z given by 

(X1,X 2 ,X 3 ,X 4 ) <~ (2/1,2/2,2/3,2/4) 

and consider the map 

given by ^ = II o ^ where II is a canonical projection. This map is surjective. In fact we have the 
correspondence 

• xi — x 3 > 0, x 2 — x 4 > 0, Xi — x 3 < x 2 — X4 

(xi,x 2 ,x 3 ,x 4 ) — (xi - x 3 ,x 2 - x 4 , 0,0) = Tp([x\ - x 3 ,x 2 - x 4 ], 0). 

• Xi — x 3 > 0, x 2 — x 4 > 0, xi — x 3 > x 2 — x 4 

(xi,x 2 ,x 3 ,x 4 ) ~ (0,0, x 3 - xi,x 4 - x 2 ) = ^(0, [x 3 - xi,x 4 - x 2 ]). 

• xi — x 3 > 0, x 2 — x 4 < 

(xi,x 2 ,x 3 ,x 4 ) ~ (0,x 2 - x 4 ,x 3 - xi,,0) = ^(0, [x 3 - xi,x 4 - x 2 ]). 

• xi — x 3 < 0, x 2 — x 4 > 

(xi,x 2 ,x 3 ,x 4 ) ~ (0,x 2 - x 4 ,x 3 - xi,,0) = ^([x 3 - xi,x 2 - x 4 ],0). 
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• x\ — x 3 < 0, x 2 — X4 < 0, x\ — x 3 > x 2 — X4 

(x 1 ,X2,x 3 ,x 4 ) ~ (xi - x 3 ,x 2 - x 4 ,0,0) = Tp(0, [x 3 - xi,x 4 - x 2 ]). 

• x\ — x 3 < 0, x 2 - x 4 < 0, x\ — x 3 < x 2 — Xi{xi,x 2 , x 3 , x 4 ) ~ (0, 0, x 3 — x\, x 4 — x 2 ) = Tp([x 3 — 
xi,x 2 - ar 4 ],0). 

This correspondence defines a map 

ip : M — >M. 

In the following, to simplify notation, we write Tp instead of Tp. 
Definition 1 For any X, X' e IR, we put 



X • X' = ip(Tp{X) •Tp{X')). 
This multiplication is distributive with respect the the addition. In fact 



(X 1 + X 2 ) • X' = ^p(lp(X 1 + X 2 ) • lp(X')). 

Suppose that Tp{X\ + X 2 ) ^ ^(Xi) + Tp(X 2 ). In this case this means that Tp(X\) + Tp(X 2 ) ^ ImTp. But by 
construction Tp(X\ + X 2 ) € ImTp and this coincides with Tp(X\ + X 2 ). For numerical application of this 
product, see ([?, GN-R] . 

2 The module gl(n,W) 

Let gl(n, ER) be the set of square matrices of order n whose elements are in IK. A matrice of gl(n 7 lR) 
is denoted by 

A = (Xa). . , 

with Xij = (Kij,0) or (0,Kij). It is clear that gl(n,IW) is a real vector space. We define a product on 
it puting 

A-B = (X ij )-(Y ij ) = (Z ij ) 

n 

with Zij — ■ Y ik . This last product being the associative product on M. Thus gl(n, IR) is an 

k=l 

associative algebra. 

Definition 1 A matrice A e gl(n,W.) is called inversible if its determinant, computed by the Cramer 
rule, is an inversible element in IR. 

Recall that the group IR of inversible elements contain 

Xi = (KuO) or {0,Ki) 

with ^ Ki. To compute the determinant, we use the classical formula of Cramer. 

Example 1. Let us consider the matrix 



M 



[1,2] [-1.3] 
[-1,3] [1,2] 



Thus 



detBi = ([1,2],0)([1,2],0)\([-1,3],0)([-1,3],0) 

= ([l,4],0)\([-3,9],0) 

- ([0,[-4,5]) 

= \([-4,5],0). 



As ([— 4, 5],0) is not an inversible element of WL, the matrix B\ is not inversible. 
Example 2. Now if 

_ / [1,2] [-1.3] \ 
B2 ~{ [-1,3] [1,7] J 

then, by the similar computation, we obtain 

detB 2 = \([-7,-4],0) 

and B 2 is invertible. 

Definition 2 If A is an invertible matrix on gl(n,WL), the inverse matrix A^ 1 of A is given by 

A-A- l = Id 

where 



Id 



( 1 ••• \ 
1 ••• 



V ••• 1 J 



with 1 = ([1,1], 0) and0= ([0,0], 0). 

The determination of A -1 can be computed using the classical rules. 
Example. If we consider the invertible matrix B 2 , we obtain 



L 7'4 J V M-1,3] [1,2] 
Let us verify that B 2 B 2 ^ 1 = Id. Using the product on M wc obtain 



B 2 B 2 -[-, -]^ hlj3] 



[1,7] M-1.3] 
\[-l,3] [1,2] 



The coefficient in place (1,1) is 



-1 1, 



o 11 = l -,-]([l,2][l,7] + [-1.3](\[-l,3])). 
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^From the definition of the product (see section 1), this element is 



an = (Ho,0)((l,2 1 0,0)(l 1 7 ) 0,0)-(0,3,1,0)(0,3,1,0) 

i 1 0,0)((1,14,0,0)-(0,10,6,0)) 

^,0,0)(l,4,-6,0) 

^,0,0)(7,4,0,0) 
1,1,0,0) 

which corresponds to [1,1]. Similarly we have a\2 = a 2 i = (0,0,0,0) and 022 
B 2 B^ =Id. 



(1,1,0,0). Thus 



3 Diagonalization 



3.1 Eigenvalues and central eigenvalues 



Let A be in gl(n, M). An eigenvalue of A is an element X e M such that there exists a vector V^Oe 
with 

A- t V = X- t V. 

Thus A" is a root of the characteristical polynomial with coefficients in the ring IR 

C A {X) = det(A - XI) = 0. 

Example. Let 

r ( [1,2] [1,2] 

[1,3] [2,5] 



We have 



and 



B *- XI -{ [1,3] [2,5] x 



X 



det(B 3 -XI) = ([1,2]\*)([2,5]\*)-[1,3][1,2] 

= [2, 10] - X[2, 5] - X[l, 2] + (\X)(\X) - [1, 6] 
= [\X){\X)-X[3,7] + [1,A]. 

Let X — ([x, y],0). It is represented in A4 by (x, y, 0, 0) or (0, y, x, 0) or (0, 0, x, y) = — (x, y, 0, 0). 

First case: det(B 3 -XI) = {x 2 , y 2 , 0, 0) - (3x, 7y, 0, 0) + (1, 4, 0, 0) = (x 2 - 3x + 1, y 2 - 7y + 4, 0, 0). 
Then dct(B 3 - XI) =0 implies 

f x 2 - 3x + 1 = 0, 
1 y 2 -7y + 4 = 0, 



G 



that is 



3± V5 
7± V33 



We obtain 

* = ([^,^!],o). 

Second case: det(S 3 - A"/) = (0, y 2 + x 2 , 2xy, 0) - (0, 7y, 7x, 0) + (1, 4, 0, 0) = (1, y 2 + x 2 - 7y + 
4, 2xy - 7x, 0). Then det(B 3 -XI) =0 implies 

1 - 2a;y + 7x = 0, 
y 2 + x 2 - 7y + 4 = 0. 

This gives 

4y 4 - 56y 3 + 261y 2 - 455y + 197 = 0. 

We have the following solutions 

(x; y) = {(-2, 8; 3, 32), (2, 9; 3, 67), (-0, 17; 0, 63), (0, 17; 6, 37)}. 
We obtain the eigenvalues 

X 4 = ([-2,8, 3.32], 0), 
X 5 = ([-0.17, 0.63], 0). 

Third case: det(B 3 -XI) = (x 2 , y 2 , 0, 0) + (Zx, 7y, 0, 0) + (1, 4, 0, 0) = (x 2 + 3x + 1, y 2 + 7y + 4, 0, 0). 
Then det(B 3 - XI) = implies 

-3± y/5 



x = 



y = 



2 ' 
-7±V33 



then 



We obtain six eigenvalues. 



^6 = ([" 



V5 



,0). 



Remark. To compute the interval-eigenvalues of a matrix A, we have to find the roots of the character- 
istical polynomial of A. But this polynomial is with coefficients in IM (or Aa) and this set is not a field 
neither a factorial ring. Then it is natural to meet some special results (e.g if we consider the second 
degree polynomial X 2 — 1 with coefficients in J| which is not factorial, it admits four roots, 1, 3, 5, 7.) In 
our example we finds 6 roots. Now if we consider the real matrix whose coefficients are the centers of 
interval-coefficients of B 3 , that is 

' 1.5 1.5 
2 3.5 
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then the eigenvalues of cb 3 are 4.5 and 0.5 which are closed to the center of X\ and X 3 . We call these 
eigenvalues, the central eigenvalues. 

Definition 3 Let A be a matrix in gl(n, WL). Lat A c be the real matrix whose elements are the center 
of the intervals of A. We say that an eigenvalue of A is a central eigenvalue if its center is (close to) 
an eigenvalue of A c . 



Remark. The determination of negative eigenvalues that is of type (0, K) is similar. Nevertheless we 
have to consider only matrices with positive entries thus we studies only the positive eigenvalues. The 
negative eigenvalues do not correspond to physical entities. 



3.2 Eigenvectors, eigenspaces 

Now we will look the problem of reduction of an interval matrix. Recall that the characteristical poly- 
nomial is with coefficient in a non factorial ring. This is the bigest change with respect the classical real 
linear algebra. 

Definition 2 Let A a square matrix with coeff dents in IM. If X is an eigenvalue of A, then every vector 
V £ M satisfying A*V = X*V is an eigenvector associated with X. 

Let Ex be the set 

E x = {V e IS" such that A*V = X f V}. 
Then Ex is a R-subspacc of M™ where n is the order of the matrix A. It is also a M submodule of IR". 

Proposition 1 Let X\ and X2 be two distinguish eigenvalues of A. Then E Xl H Ex 2 = {0}. 

Proof. Let V be in E Xl n Ex 2 - We have 

A l V = X t V, 
A*V = X 2 V. 

This X X V \ X 2 V = (X 1 \ X 2 )V = 0. As II is without zero divisor, we have X x \ X 2 = or V = 0. We 
deduce E Xl fl Ex 2 = {0}. 



Proposition 2 Let Ca{X) be the characteristical polynomial of A. If the real polynomial C c A {X) asso- 
ciated with the central matrix of A is a product of factor of degree 1, then Ca(X) admits a factorization 
on M 

We have seen that Ca(X) can be have more than degree^A^) roots. If Xi, ■ ■ ■ ,X n are the central 
roots, we have the decomposition 



C A {X) =a n ]J(X\X i ). 



i=i 

Example. If we consider the matrix 



[1.2] [1,2] 
[1,3] [2,5] 
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then Cb 3 (X) admits X\, ■ ■ ■ ,X 6 as positive roots. The central eigenvalues are X\ and X 3 and we have 

det(B 3 \ XI) = (X \ X X ){X \ X 3 ). 

If we consider the roots X 2 = ([ — - — , ],0),and if we assume that C'b 3 s( X) = (X \X 2 )(X \Y), 



we obtain 



Y = (3,7, 



2 ' 2 

3-\/5 7 + V33 N ,3 + V5 7 - V33 



2 ' 2 

which does not correspond to a positive eigenvalue. 



) = (- 



,0,0) 



Theorem 1 For any n-ule of roots (X\, ■ ■ ■ ,X n ) such that Ca(X) = a n J| (X \ Xj), and if for 

i=l 

i = ,n the dimension of Ex i coincides with the multiplicity of Xi, then we have the vectorial 

decomposition M = Q>i e iEx i where the roots X i7 i € / are pairwise distinguish. 



Example. Let us compute the eigenspaces of B 3 associated to the central eigenvalues. 

v i ^ + ^> 7+V33, 
• x i = ([ — ^ — . 7i J' )- 



Let V = 



2 ' 2 

v 2 



£ IE. Then 



(A - X X I) 



V 2 



is equivalent to 



2 ' 2 
V ([l-3],0)Vi 



(0, 



-1 + V5 -3 + ^ 



= 



w ) 



that is 



This gives 



M l±^,i±^]v I + [1 , 21 v 2 = o, 



2 2 

r l + V5 3 + V33 



K 2 



[1.2] 



If we choose Vi = ([1, 1], 0) we have 



V 2 = 



1 + V5 3 + V33 

[1.2] 

-1 + V5 -3 + V33 



]-N-i,^]) 



x([ _3 ± ^33 ) _l^I ]) 
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Thus the -eigenvectors are of the form 

V = 



([1,1], 0) 



x([ _M^,_i±^],o) 



Remark. We can choose V\ such that all tje coordinatez of V are positive. For example if V\ = [1,2] 

([1,2],0) 
then V = | 1 + V5 3 + ^ 

([ 2 , 2 J.0) 



Let V 



Y) I e M. Then 

V2 



(^-^)( £ )=o 



is equivalent to 



(F^^^Wr ([l,2],0)V a 



([1-3], 0)^ 



([ 



0)V 2 J 



that is 



This gives 



V 2 - 

If we choose Vi = ([1, 1], 0) we have 



■([- 



2 ' 2 
-1 + V5 -3 + ^ 



])^i 



V 2 



[1.2] 

-1 + V5 -3 + ^33 
2 ' 2 1 

[1,2] 

-1 + V5 -3 + V33 n 



W-i,^D 



Thus the A^-eigenvectors are of the form 

V = 



([1,1], 0) 
([^,^],0) 
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4 The Exponential map 

We define the exponential map 

Exp : gl(n,W) — > gl(n,W) 
in a classical way by series expansions. If the matrix A is diagonalizable, then 

D = P- 1 AP 

is diagonal and Exp(A) is a diagonal matrix whose diagonal element are the exponential of the eigen- 
values. 

Example. Let 



The central eigenvalues are 



„ _ I [1.2] [1,2] 
{ [1.3] [2,5] 



*-([^,^l,0l 
* = ([^.^],0). 



and we have 



with 



We deduce 



D 



([1,1], 0) 



X x 
X 3 



([1,1], 0) 



Exp(B 3 ) = P. 



{[exp( 3 + 2 ^ ), exp{ 7 + 2 V ^ ], 0)) 









{[exp{ 3 ^ ),exp( 7 2 V ^ ],0)) 



In a forthcomming paper, we apply this calculus to solve linear differential system. 



References 



[1] Goze Nicolas, Remm Elisabeth. An algebraic approach to the set of intervals. arXiv math 0809.5150 
(2008) 



11 



